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We present results of systematic fully relativistic first-principles calculations of the uniaxial mag- 
netic anisotropy energy (MAE) of a disordered and partially ordered tetragonal Fe-Co alloy using 
the coherent potential approximation (CPA). This alloy has recently become a promising system for 
thin ferromagnetic films with a perpendicular magnetic anisotropy. We find that existing theoretical 
approaches to homogeneous random bulk Fe-Co alloys, based on a simple virtual crystal approxima- 
tion (VCA), overestimate the maximum MAE values obtained in the CPA by a factor of four. This 
pronounced difference is ascribed to the strong disorder in the minority spin channel of real alloys, 
which is neglected in the VCA and which leads to a broadening of the d-like eigenstates at the Fermi 
energy and to the reduction of the MAE. The ordered Fe-Co alloys with a maximum Llo-like atomic 
long-range order can exhibit high values of the MAE, which, however, get dramatically reduced by 
small perturbations of the perfect order. 

PACS numbers: 71.15.Rf, 71.23.-k, 75.30.Gw, 75.50.Bb 



I. INTRODUCTION 

The binary ferromagnetic Fe-Co alloy has been known 
for a long time as a system with the maximum 
spontaneous magnetization among the transition-metal 
systemsii Its basic magnetic properties, such as the con- 
centration dependence of the alloy magnetization (the 
Slater-Pauling curve) in the ground-state body-centered 
cubic (bcc) structure, have been reproduced successfully 
by ah initio electronic structure calculations in a num- 
ber of studies.^""* This system attracted renewed interest 
several years ago, after a theoretical prediction of a giant 
uniaxial magnetic anisotropy energy (MAE) of the body- 
centered tetragonal (bet) Yei-xCox alloysi^ The uniaxial 
MAE is defined quantitatively as the difference of total 
energies for the magnetization direction parallel to the 
tetragonal a and c axis, = E^'^°^'^ - E^^°'^\ The cal- 
culated MAE reached high values, w 800 /ieV/atom, 
obtained however only in a narrow range of the Co con- 
centration, 0.55 < a; < 0.65, and of the tetragonal strain, 
1.2 < c/a < 1.25. The combination of the high MAE 
with the high magnetization makes the bet Fe-Co system 
a promising material for fabrication of ferromagnetic thin 
films with a perpendicular magnetic anisotropy (the mag- 
netic easy axis perpendicular to the film plane), which 
might be relevant as high-density magnetic recording me- 
dia. 

Most of experimental realizations of this tetrago- 
nal system employed the possibility to grow epitaxially 
strained Fe-Co alloy films of a varying composition on 



different non-magnetic transition-metal substrates with 
the face-centered cubic (fee) structure, such as Pd(OOl), 
Rh(OOl), and Ir(OOl),^ The different lattice parame- 
ters of the substrates enable one to scan a relatively 
wide interval of the c/a-ratio of the Fe-Co films, namely 
1.13 < c/a < 1.24 (note that c/a = 1 and c/a = \/2 
refer respectively to the ideal bcc and fee lattices). The 
measurements of the in situ magneto-optical Kerr effect^ 
(MOKE) confirmed the theoretical prediction of Ref. 5 on 
a qualitative level; in particular, a strong perpendicular 
magnetic anisotropy of the films was observed for compo- 
sitions and tetragonal distortions in a rough agreement 
with the calculated trends of the bulk K^. However, sys- 
tematic quantitative experimental studies of the MAE of 
the films have not been performed yet; available values 
of the for a few selected systemsSrJS are appreciably 
smaller than the calculated ones. The missing informa- 
tion on the MAE has partly been compensated by mea- 
surements of X-ray magnetic circular dichroism (XMCD) 
spectra, which provide orbital magnetic moments of both 
alloy constituents.^ The latter measurements, performed 
on the Rh(OOl) substrate corresponding to c/a w 1.24, 
revealed enhanced orbital moments in Fei-^jCoa; films 
around x ~ 0.6, in a close relation to the predicted max- 
imum of the bulk uniaxial MAE. 

Nevertheless, the existing agreement between the the- 
ory and the experiment should be taken with caution. 
The original approach of Ref. ^ employed the so-called 
virtual crystal approximation (VCA), in which the true 
species of a random binary alloy are replaced by a sin- 
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gle element with an effective atomic number given by the 
alloy composition. The VGA has been implemented in 
various ab initio methods and used for a number of al- 
loys of neighboring elements in the periodic tablej^i^iii^ 
Its application to the bcc Fe-Co and fee Co-Ni systems^ii^ 
yields concentration trends of the spin magnetic moment 
in a very good agreement with experiment;^ the orbital 
magnetic moments seem to be described reliably in the 
VGA as well. However, a recent first-principles study 
of the random bet Fei_a;Coa; alloys, based on a super- 
cell technique applied to a few special Co concentrations 
(x = 0.5, 0.625, and 0.75), has shown that a realistic 
treatment of the chemical disorder can reduce the MAE 
by a factor of 1.5—3 as compared to the simple VGA. 

The main purpose of this work is to present and dis- 
cuss the results of a systematic ab initio calculation of the 
MAE for the disordered bulk tetragonal Fe-Co systems 
obtained by using the coherent potential approximation 
(CPA) as a basic tool of the theory of metallic alloys . ^^'"'^^ 
For homogeneous bet alloys, we compare results of the 
VGA and the CPA, show the big difference between them, 
and identify the underlying physical mechanism in terms 
of the electronic structure. Moreover, we investigate the 
effect of a complete and an incomplete atomic long-range 
order on the MAE; this study is motivated by the exist- 
ing prediction of a large uniaxial MAE in stoichiometric 
perfectly ordered tetragonal FeCo and FeGos systems. 
Similar theoretical studies have so far been done mainly 
for tetragonal FePt alloysii"— and very recently also for 
the FeCo alloy,— These works are confined to stoichio- 
metric equiconcentration alloys; in general, their results 
reveal a decrease of the MAE due to imperfect chemical 
ordering. 



II. MODELS AND COMPUTATIONAL DETAILS 

All calculations of this study employed a bet lattice 
with the Wigner-Seitz s radius equal to that of pure bcc 
iron, s = 2.662 aq where oq = 5.292 x 10~^^ m denotes 
the Bohr radius. The neglect of volume relaxations on 
the MAE in a broad range of lattice parameters and of al- 
loy concentrations is justified by a more general study of 
Ref. [13 within the VGA. This structure model — the vol- 
ume conserving tetragonal distortion (the Bain path) — 
coincides with that used in the original study"" and in the 
recent study of the partially ordered FePt alloy .J ^ 

The effect of atomic long-range ordering has been stud- 
ied for the case of LIq (CuAu) order relevant for com- 
positions not far from the equiconcentration Fco.sCoo.s 
system. The bet structure is partitioned in two simple 
tetragonal sublattices (alternating atomic planes perpen- 
dicular to the tetragonal c axis). For the Fei_2;Goa; sys- 
tem, an additional concentration variable y is introduced, 
such that < y < min{a;, 1 — a;}, which defines the chemi- 
cal composition of both sublattices: Fei-x+yCox-y in the 
Fe-enriched planes and Fei-x-yCox+y in the Co-enriched 
planes. Note that the homogeneous solid solution corre- 



sponds to y = 0, while the other end of the y- interval 
describes the maximum LIq order compatible with the 
given total Go concentration x. This model is a natural 
generalization of the model used for the stoichiometric 
FePt systems Ji"— 

The eletronic structure calculations were done by 
means of the tight-binding linear muffin-tin orbital (TB- 
LMTO) method in the atomic sphere approximation 
(ASA)^i^ with a full inclusion of relativistic effects by 
solving the one-electron Dirac equation^'^'^^ with the spd- 
basis of the valence orbitals. The effective spin-polarized 
potential was constructed in the local spin-density ap- 
proximation (LSDA) using the parametrization of the 
exchange-correlation term according to Ref. 25. The 
effects of alloying were treated in the CPA for all 
systems;'22'24 foj- the homogeneous random alloys (y = 0), 
the simple VGA was used as well. 

The reliable evaluation of the MAE in metallic systems 
represents a difficult task, especially for cubic 3d transi- 
tion metals owing to the weakness of the spin-orbit in- 
teraction and the high symmetry of cubic structures i^>^ 
The situation is more favorable for uniaxial systems 
(tetragonal, trigonal, hexagonal) ^^'^^ and for systems 
with /-electrons, see Ref. |30| for a review. In a number of 
existing studies, including the very recent onesr' i^^'^^ the 
magnetic force theorem'*^ is used and the total energy dif- 
ference K^i is approximated by the change in the sum of 
occupied valence one-particle eigenvalues. In this study, 
we go beyond this approximation and evaluate the MAE 
directly from the total energies of the fully self-consistent 
solutions for both magnetization directions which 
requires high accuracy in total-energy calculations. We 
have used uniform sampling meshes of about 10"' k-points 
for averages over the full Brillouin zone (BZ) of the bet 
lattices, while about 5 x 10"* k-points have been used for 
the BZ of the simple tetragonal Bravais lattices of the 
LIq ordered systems. The total energies were converged 
to 0.1 /ieV/atom for all systems; this accuracy is suffi- 
cient in view of the resulting values of (see below). 

III. RESULTS AND DISCUSSION 
A. Homogeneous random alloys 

1. Magnetic anisotropy energy 

The dependence of the uniaxial MAE of the random 
Fei-^^Co^: alloys on the chemical composition and the 
c/a-ratio, calculated in the VGA and in the CPA, is 
shown in Fig. [1] for better transparency, only the pos- 
itive values of the are displayed in the plots, while 
the cases with negative Ku have been omitted (marked by 
white color). The VGA results (Fig.[T^) agree nicely with 
both previous studiesf^*^ in particular, a sharp maxi- 
mum of the MAE, « 809 /xeV/atom, obtained for 
X — 0.6 and c/a = 1.24, is clearly visible. The maxi- 
mum MAE of Ku ~ 800 /xeV/atom for x = 0.6 and for 
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c/a between 1.22 and 1.24 was reported in Ref . li' (with a 
slight sensitivity to the particular model of the volume re- 
laxation employed) and practically the same values were 
obtained in the original VGA study.^ Minor quantitative 
differences might be ascribed to different computational 
schemes employed. The MAE in the CPA (Fig. [Hd) ex- 
hibits a similar trend as in the VGA; however, the CPA 
maximum is significantly smaller, iCu ~ 183 /leV/atom. 
The latter value is obtained for x = 0.6 and c/a = 1.23, 
but the maximum is very flat and essentially the same 
values can be found in a broad region of 0.5 < x < 0.65 
and 1.21 < c/a < 1.28, see Fig. [T]d. The maximum 
in the CPA is reduced by a factor of four as compared 
to that in the VGA, which is an even stronger reduction 
than that reported in Ref. Q 

The small MAE values in the CPA might be used to ex- 
plain the difference between the large values predicted 
in the VGA and the much smaller values inferred from 
measurements on thin Fe-Go films.®"— Undoubtedly, the 
better description of the chemical disorder in the CPA 
as compared to the VGA contributes to the reduction of 
the MAE of prepared alloy thin films. However, a thor- 
ough quantitative analysis of this point is impossible at 
present and it could even be misleading because of the 
well-known uncertainty of the LSDA to yield quantita- 
tively correct MAEs in 3d transition-metal systems. For 
this reason, an analysis of the trends seems to be more 
appropriate, see, e.g., Ref. [H and references therein. 

In general, the LSDA underestimates both the MAEs 
and the orbital magnetismJ^iiii^!^ A counterexample 
to this rule is represented by the MAE of the ordered 
FePt alloyr^^ which is overestimated in the LSDA prob- 
ably owing to an interplay of intraatomic Coulomb cor- 
relations, weak exchange splitting and strong spin-orbit 
interaction of Pt atoms. The present results for the bet 
Fei_a;Coa: alloys seem to confirm the general rule as can 
be documented by two facts. First, the effective MAE 
of strained Fe-Go thin films, given by the difference of 
the bulk and the magnetostatic shape anisotropy en- 
ergy Ksh, comes out positive for alloys in the vicinity of 
X 0.6 and c/a 1.25 (near the maximum of in the 
CPA) . This situation would thus lead to an out-of-plane 
orientation of the thin-film magnetization, in a qualita- 
tive agreement with existing experimentsii"— Note that 
Ksh — (a^o/2)M^/F, where fio is the permeability of the 
vacuum and where M and V denote, respectively, the al- 
loy magnetic moment per atom and the atomic volume. 
In the present case, M w 2.14 fiB, where /ie is the Bohr 
magneton, one obtains Kgh ~ 132 /^eV/atom, which lies 
below the bulk MAE (K^ w 180 ^eV/atom). However, 
for the case of x — 0.5 and c/a = 1.13, corresponding 
roughly to the Feo.sCoo.s films grown on the Pd(OOl) 
substrate, our calculated values are — 124 //eV/atom, 
M « 2.27 /iB, and Ksh ~ 148 /ieV/atom, which indi- 
cates an in-plane orientation of the magnetization, in 
contrast to the out-of-plane orientation observed at low 
temperatures f2, Second, the calculated orbital magnetic 
moments in the bulk Fei_3;Go2; alloys, plotted in Fig. [5] 
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FIG. 1. (Color online) The uniaxial MAE of random bet 
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in the CPA (b). Only positive values of the Ku are dis- 
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FIG. 2. (Color online) The orbital magnetic moments in the 
random bet Feo.4Coo.6 alloy with magnetization along z axis 
as functions of the c/a-ratio: the alloy orbital moment in the 
VCA and the species-resolved orbital moments in the CPA. 
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for X = 0.6, are appreciably smaller than the values ob- 
tained from the measured XMCD spectra by employ- 
ing the sum rulesj^ In particular, the measured val- 
ues for films with x = 0.6 and c/a = 1.24, namely, 
M^-i, = 0.21 ± 0.03 //B and M^" = 0.32 ± 0.04 ^b, exceed 
their bulk theoretical counterparts by a factor of three. A 
similar relation of the experiment and the LSDA theory 
(quantified roughly by a factor between 1.5 and 2) was 
obtained for bcc Fe-Co alloysi^ii^'^ Moreover, the data 
in Fig. [2] prove that the VGA and the CPA yield quali- 
tatively different trends of the orbital magnetism versus 
the c/a-ratio: the pronounced maximum at c/a = 1.24 
in the VGA is replaced by very flat maxima of both lo- 
cal orbital moments in the GPA. This difference indicates 
that applicability of the VGA to all details of magnetism 
of the Fe-Go system is limited, being confined probably 
only to the cubic structures. '^'^'^ 

In view of the above mentioned problems to evaluate 
reliably the MAEs and the orbital magnetic moments, we 
have not attempted to recalculate our results by tech- 
niques going beyond the LSDA^^i^^i^i^ but have fo- 
cused on the strong difference between the VGA and 
the GPA and on the role of disorder on the MAE in the 
tetragonal Fe-Go systems. 
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FIG. 3. (Color online) The spin-polarized local densities of 
states of the random bet Feo.4Coo.6 alloy with c/a = 1.25 as 
functions of the energy. The dotted vertical line denotes the 
position of the Fermi level. 



2. Electronic structure 

The giant MAE obtained in the VGA was ascribed 
originally to the changes in the band structure accom- 
panying the tetragonal distortion of the effective FeGo 
crystal.— The explanation rests on a coincidence of two 
particular eigenvalues at the F point, which occurs in 
the minority spin (spin-down) channel. For a special al- 
loy composition and a special c/a-ratio, this coincidence 
takes place just at the Fermi energy corresponding to the 
particular filling of the spin-down band. The large con- 
tribution of this eigenvalue pair to the MAE can be un- 
derstood in terms of the second-order perturbation the- 
ory, in which the effect of the weak spin-orbit interaction 
can be safely included. In this approach, the enhanced 
value of the of the bulk Fe-Go alloys can be explained 
due to the coupling between d^^-yi and dxy states me- 
diated by the orbital momentum operator L^, see Ref. 
for details. The physical origin of the perpendicular mag- 
netic anisotropy of the strained Fe-Go films is thus simi- 
lar to that of Au covered Go monolayers on an Au(lll) 
substratei^ In the latter case, however, the resulting pos- 
itive MAE is caused by the coupling between dxz and 
dyz states. 

The usual band structures are not relevant for the 
electronic structure of random alloys, for which more 
appropriate quantities, such as the densities of states 
(DOS) and the Bloch spectral functions (BSF), have to 
be studied.— li^i^ Motivated by the above explanation of 
the giant MAE, we present these quantities for the ran- 
dom bet Fei-xCOx alloys in the GPA, evaluated without 
the spin-orbit interaction, i.e., in the scalar-relativistic 



approximation . 

The spin-polarized local DOSs are shown in Fig. |3]for 
the system with x = 0.6 and c/a — 1.25. The shapes 
of the individual DOSs prove a very weak disorder in 
the majority spin (spin-up) channel, whereas the spin- 
down band exhibits a regime of a stronger scattering. 
This type of spin-dependent disorder was found in the 
bcc Fe-Go system a long time ago^^ it is responsible, 
e.g., for the observed concentration trend of the residual 
resistivity.^ Note that the strong disorder is present in 
the spin-down d-band, which is only partially occupied 
and which thus contributes a lot to the alloy total energy 
and, consequently, to the MAE. 

The spin-polarized BSFs at the F point for the 
Feo.4Goo.6 alloy with three tetragonal distortions are 
plotted in Fig. The BSFs were resolved according to 
the irreducible representations of the point group 
of the bet lattice.— The odd (ungerade) representations, 
belonging to p orbitals, have been omitted in the plots 
since their contribution is negligible in the displayed en- 
ergy interval. The relevant even (gerade) representations 
are: Aig (orbitals s and d^-i), Big (dx-^-y-^), (dxy), 
and Eg (dxz and dyz). One can see that all spin- up BSF's, 
consisting of narrow Lorentzian peaks, can be interpreted 
as Bloch-like eigenstates with finite lifetimes due to a 
weak disorder. The spin-down BSF's above the Fermi 
energy (mainly ^ig-like) reflect effects of strong disorder 
(non-quasiparticle behavior); the profiles at and below 
the Fermi energy are Lorentzian peaks again. However, 
their widths are clearly bigger than the spin-up widths, in 
agreement with the spin-dependent disorder manifested 
in the DOS. 
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FIG. 4. (Color online) The spin-polarized symmetry-resolved 
Bloch spectral functions of the random bet Feo.4Coo.6 alloy 
at k = r as functions of the energy for three values of the 
c/a-ratio. The dotted vertical lines denote positions of the 
Fermi levels. For details, see the text. 



The maximum of the MAE in the VGA is found for an 
alloy with a: = 0.6 and with the tetragonal strain close 
to c/a = 1.25; its spin-down BSF in the CPA (Fig. H 
middle panel) shows an analogy of the coincidence of the 
Sig (d^i_yi) and iAxy) levels at the Fermi energy. 
However, the disorder-induced smearing of both peaks is 
at least as big as their separation, which suppresses the 
contribution of this eigenvalue pair to the MAE in the 
framework of the second-order perturbation theory. This 
feature proves that the minority-spin Fe-Co disorder is 
strong enough to reduce significantly the MAE, which 
explains the very big difference between the VGA and 
the CPA results, see Section HTl A II 



B. Partially ordered alloys 

The adopted model of the partially ordered tetrag- 
onal Fei_a;Coa; alloys, namely, the model of two sub- 
lattices with chemical compositions Fei_2;+yCoa;_y and 
Fei_j;_yCoa;+y (Section[n| and the use of the CPA (in 
contrast to supercell techniques) enable one to study the 
MAE as a function of three continuous variables: x, y, 
and c/a. Such a full three-dimensional scan is, however, 
computationally very demanding; we have thus confined 
our study to Co concentrations close Xox — 0.6, for which 
the maximum MAE was found in the homogeneous bet 
alloys both in the VCA and in the CPA. The depen- 
dence of the calculated on the tetragonal distortion 
for X — 0.6 is shown in Fig. [5] for three degrees of the 
LIq order. One can see that the very flat low maximum 
iiTu ~ 183/ieV/atom for the completely disordered alloy 
(y — 0) treated in the CPA is replaced by a significantly 
enhanced pronounced maximum « 450 /xeV/ atom ob- 
tained for the alloy with the maximum order (y — 0.4) 
compatible with the given Co concentration {x — 0.6). 
The latter maximum is about two times smaller than 
that in the VCA {K^ « 809 /ieV/atom), which can be 
ascribed to the fact that the system with x = 0.6 and 
y — 0.4 has one sublattice without disorder (pure Co), 
whereas the other one is disordered, containing the rest 
Co and all Fe atoms. The MAE of the system with an in- 
termediate degree of the Llg order (y = 0.2) exhibits the 
maximum and the trend very close to the fully disordered 
alloy. 

The detailed dependence of the on the concentra- 
tion variable y is presented in Fig.|6]for two cases: the sto- 
ichiometric alloy {x — 0.5, c/a ^ 1.24) and the previous 
off-stoichiometric alloy {x — 0.6, c/a — 1.27). The par- 
ticular values of the c/a-ratio chosen in both cases corre- 
spond to the maximum value obtained for alloys with 
the maximum order, i.e., for x = 0.6, y = 0.4 (Fig. [5]) and 
for X = y = 0.5. One can see that the perfectly ordered 
{y — 0.5) stoichiometric alloy leads to a very high MAE of 
Ku ~ 580 /j,eV/atom. This value agrees well with Ref . [T^. 
where the maximum MAE for the same LIq ordered sys- 
tem was obtained as ~ 520/ieV/atom. The convex 
shape of both dependences in Fig. [S] proves that the very 
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FIG. 5. (Color online) The uniaxial MAE of the bet 
Feo.4Coo.6 alloy as a function of the tetragonal strain c/a for 
different degrees of the Llo-like atomic order: the maximum 
order (j/ = 0.4, diamonds), an intermediate order (j/ = 0.2, 
full circles), and the complete randomness (y = 0, crosses), 
all treated in the CPA. For a comparison, the MAE of the 
completely random alloy in the VGA is displayed as well (tri- 
angles) . 



high MAEs can be obtained only for systems with the 
maximum order; even a small amount of additional dis- 
order is detrimental to the MAEs and reduces them to 
much lower values of the completely random alloys. ^'^ 

Strong sensitivity of the MAE to the degree of the LIq 
atomic order has recently been reported for the stoicho- 
metric FePt alloy on an fee lattice with small tetrago- 
nal distortions The results of Ref. [l^ indicate that the 
chemical ordering is a more important factor for high 
MAE values than the tetragonal distortion. Our results 
for the Fe-Co system witness that both factors are of 
equal importance, see Fig. [Sj We believe that validity of 
this type of conclusions depends also on the range of rele- 
vant variables: the tetragonal distortion was varied over a 
narrower interval (0.94 < c/a < 1.06) in Ref. [l9l. whereas 
a wider interval (1.15 < c/a < 1.35) has been covered 
in our study. Another difference between the two sys- 
tems lies in the dependence of the MAE on the c/a-ratio: 
the MAE is an ever increasing function of c/a in FePt)^ 
in contrast to the maxima found for the Fe-Co systems 
around c/a w 1.25, see Fig. [T]and Fig. [S] These trends 
reflect probably the different origins of the high MAE 
in these systems. In the FePt alloys, the iron sites are 
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FIG. 6. The uniaxial MAE of two bet Fei-^^Coj, alloys as a 
function of the degree of the Llo-like atomic order: for x = 0.5 
and c/a — 1.24 (full circles), and for x — 0.6 and c/a — 1.27 
(open squares) . 



responsible for strong exchange fields and the platinum 
sites provide strong spin-orbit interaction, both effects 
being only little dependent on the tetragonal distortion 
of the lattice. In the Fe-Co alloys, both species are fea- 
tured by strong exchange splittings and weak spin-orbit 
couplings; the high MAEs are heavily based on collective 
properties (band structure) of the tetragonal systems. 



IV. CONCLUSIONS 

We have shown by means of first-principles LSDA cal- 
culations that chemical disorder has a strong influence 
on the uniaxial MAE of the bulk tetragonal Fe-Co sys- 
tems. First, the complete neglect of the disorder, in- 
herent to the simple VCA, overestimates the MAE by a 
large factor, whereas the more sophisticated CPA leads 
to a drastic reduction of the MAE. The latter low MAEs 
are quite close to the magnetic shape anisotropy energy 
and the resulting estimated stability of the perpendicular 
magnetic anisotropy of thin strained Fe-Co films comes 
out probably smaller than the measured one. Second, 
the disorder-induced reduction of the MAE is due to 
the strong scattering regime in the minority-spin chan- 
nel, which smears the Bloch-like eigenstates at the Fermi 
energy. Third, an analysis of the calculated orbital mag- 
netic moments points to a non-negligible underestima- 
tion of the orbital magnetism and, most probably, of the 
MAEs by the LSDA. Fourth, the Llo-hke atomic long- 
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range order leads to an enhancement of the MAE values. 
The most pronounced enhancement is obtained for the 
maximum degree of the order, while imperfect ordering 
reduces the MAE very rapidly as compared to the per- 
fectly ordered systems. Similar correlations between the 
atomic order and the MAE have recently been found in 
the FePt alloyi^"— and observed in artificially synthe- 
sized FeNi films4i We believe that this interplay should 



be taken into account in a future search of advanced ma- 
terials for high-density magnetic recording. 
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